clc;
syms u v;
syms a b c positive;
x=a*sin(u)*cos(v);
y=b*sin(u)*sin(v);
z=c*cos(u);

I=surf_integral([0,0,x*y+z],[x,y,z],[u,v],[0,pi/2],[0,2*pi])